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Few-body systems with large scattering length a have universal properties that do not depend on 
the details of their interactions at short distances. The rate constant for three-body recombination 
of bosonic atoms of mass m into a shallow dimer scales as ha 4 /m times a log-periodic function of the 
scattering length. We calculate the leading and subleading corrections to the rate constant which 
are due to the effective range of the atoms and study the correlation between the rate constant and 
the atom-dimer scattering length. Our results are applied to 4 He atoms as a test case. 
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A. Introduction 

The properties of low-energy bosons are dominated by 
S-wave interactions. If the interaction is short ranged, 
these properties can be parametrized by the effective 
range expansion. The typical range of the interaction 
also defines a natural low-energy length scale I For 
a finite range potential I is simply given by the range of 
the potential. If the potential has a van der Waals tail 
— Ce/r 6 , it is determined by the van der Waals length 
Ivdw = {mCe/h 2 ) 1 / 4 , where m is the mass of the bosons. 
For a generic physical system, the effective range param- 
eters, such as the scattering length a and the effective 
range r e , are of the same order of magnitude as I. In 
some systems, however, the scattering length a is much 
larger than I while the effective range r e ~ I is still of 
the same order of magnitude. This situation requires a 
fine tuning of the parameters in the underlying potential. 
Typical examples are 4 He atoms and nucleons where this 
fine tuning is accidental, or alkali atoms near a Feshbach 
resonance where it can be controlled experimentally by 
changing an external magnetic field. 

Physical systems with large \a\ display a number of 
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interesting effects and universal properties that are in- 
dependent of the details of the interaction at short dis- 
tances of order I The simplest example for pos- 
itive a is the existence of a shallow two-body bound 
state (called a dimer) with universal binding energy 
B 2 = Ti 2 /(ma 2 )[l + 0(l/a)]. Similar relations hold for 
other observables, which can generally be described at 
low energies by a controlled expansion in l/\a\. Par- 
ticularly interesting is the structure of the three-boson 
system. It has universal properties that include a ge- 
ometric spectrum of three-body bound states (so-called 
Efimov trimers), log-periodic dependence of three-body 
observables on the scattering length, and a discrete scal- 
ing symmetry @, H, H| . These properties can be studied 
using effective theories which provide a powerful frame- 
work to exploit the separation of scales between a and I 

QUI. 

Since we are mainly interested in applications to 
atomic systems, we will refer to the bosons as atoms in 
the following. In particular, we focus on three-body re- 
combination, the process in which three atoms collide to 
form a diatomic molecule and an atom. The energy re- 
leased by the binding energy of the molecule goes into the 
kinetic energies of the molecule and the recoiling atom. If 
the momenta of the incoming atoms are sufficiently small 
compared to 1/a, the momentum dependence of the re- 
combination rate can be neglected. The recombination 
event rate constant a is defined such that the number 
of recombination events per time and per volume in a 



2 



gas of cold atoms with number density ua is an A . If 
the atom and the dimer produced by the recombination 
process have large enough kinetic energies to escape from 
the system, the rate of decrease in the number density of 
atoms is 



—n A = -3an A 
at 



(1) 



For large positive a, recombination can go into the shal- 
low dimer with binding energy B2 = h 2 /ma 2 and an- 
other atom. If no deep dimers are present, this is the 
only channel for recombination to occur. Three-body re- 
combination into deep dimers can also be treated using 
similar methods [l|, 0, d, Q. However, here we assume 
that no deep dimers are present. We focus on a > and 
consider 3-body recombination into the shallow dimer. 



B. Theoretical framework 

We first briefly review the universal expressions at 
leading order in the expansion in I /a. Dimensional anal- 
ysis together with the discrete scaling symmetry implies 
that a is proportional to ha /m with a coefficient that 
is a log-periodic function of a with period tt/sq. To very 
high accuracy, the three-body recombination constant a 



can be expressed as 



arec-body 



ha 4 



a = 67.12 sin>o m(a/a*) + 1.67) [1 + 0{l/a)] , (2) 

m 

where sq = 1.00624. . . is a transcendental number and 
a* is the scattering length at which there is an Efimov 
trimer at the atom-dimer threshold. Note also that an 
analytic expression for the coefficient 67.12 has recently 
been derived [H, [HI]. The most remarkable feature of 
the expression in Eq. ((2]) is that the coefficient of ha A jm 
oscillates between zero and 67.12 as a function of ln(a). 
In particular, a has zeroes at values of a that differ by 
e 7r/so ~ 22.7. The locations of these zeros can be ex- 
pressed as a k (e 7r / s °)"0.20 a* where n is an integer. 
The maxima of a/a 4 occur at a p=> (e 7r / s °)"0.94 a*. These 
maxima are close to the values a — (e 7r / s °)"a >t for which 
there is an Efimov trimer at the atom-dimer threshold. In 
the latter case, atom-dimer scattering becomes resonant 
and the atom-dimer scattering length diverges. 

At leading order in I /a, the atom-dimer scattering 
length can be expressed as [H, 3 : 

a ad = (2.15cot(s ln(a/a»)) + 1.46) a [1 + 0(l/a)] . (3) 

Solving Eq. ([3]) for a*, we can eliminate a* from Eq. ^ 
and obtain an analytical expression for the universal cor- 
relation between a and aw. In this paper, we study the 
correlation between a and a ac i up to second order in in 
the expansion in I /a. To this order, these corrections 
are fully dermined by the S-wave effective range r e [l4| . 
Therefore, we refer to them as effective range corrections. 



The quantities a and a a d can both be calculated from 
the following integral equation: 1 



167TT 4 f A 

A s (p,k;E) = - 1 - V (p 7 k;E) + - / dqV(p,q;E) 

1 l r e 7T Jo 

q 2 A S (q,k;E) 

-7 + (h 2 -E- it) 1 ' 2 ~ ^ e (h 2 -E-^Y 



(4) 



where A is an ultraviolet cutoff, r e is the effective range, 
and 7 denotes the position of the bound state pole in the 
atom-atom Green's function: 



7 - (l - V 7 ! - 2r e /a) ^, 



(5) 



where we have set h = m 
potential V(p, k; E) is 



1 for convenience. The 



V(p,k;E) 



1 , fp 2 +pk + k 2 -E 
In 



2pk \p 2 — pk + fc 2 — E 



G 3 (A), (6) 



where G3 (A) represents the contribution from the contact 
three-body interaction. For r e = 0, it is given by 



G 3 (A) re -J- cot(s ln(Aa,) + 2.45) . 
A- 2 



(7) 



but for finite effective range its form is more complicated. 
In the following, we will directly calculate the correlation 
between a and a a d, therefore a» will not explicitly ap- 
pear. In our case of positive a, j 2 is the binding energy 
of the shallow dimer. At leading order (LO) in I /a, 7 
reduces to 1/a but the two quantities differ if the effec- 
tive range is included. We have chosen 7 and r e as our 
2-body inputs instead of a and r e , because this choice 
keeps the location of the dimer pole fixed which leads 
to a faster convergence of the effective-range expansion. 
The effective theory expansion is then in powers of 7L 

From the general solution of Eq. ^ one can directly 
obtain a and a ac i- The three-body recombination con- 
stant a is given by 



o 



(1 - 7^e) 



0) 



7 2 V3 (1 - jr e /2y 
while the atom-dimer scattering length is given by 



a ad = - T - As (0,0; - 7 2 ). 

6TT 



(8) 



(9) 



In principle one can obtain the effective range correc- 
tions to all orders by solving the integral equation in 
Eq. (J4j) directly. A potential problem comes from the 
dimer propagator in the second line of Eq. (H|) . In addi- 
tion to the pole from the shallow dimer with 7 ~ 1/a, it 



1 For a derivation of this equation from effective field theory, see 
e.g. Ref. Q. 
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also has a deep pole with 7 ~ l/r e . This pole is an arte- 
fact of the form of the dimer propagator and is outside 
the range of validity of the effective theory. For nega- 
tive effective range, the pole is on the unphysical sheet 
of the complex momentum plane and causes no problems 
in solving Eq. (j4]). For positive effective range, the pole 
is on the physical sheet and leads to instabilities in the 
three-body equations for cutoffs of the order ~ l/r e and 
larger. We avoid this problem by expanding the kernel 
of Eq. ([3]) in the effective range r e and solving the result- 
ing integral equation. This allows to calculate the range 
corrections perturbatively up to order {"fl) 2 . For higher 
orders, things become more complicated and a new three- 
body parameter enters [H]]. Moreover, for particles with 
an interaction displaying a van der Waals tail ~ Ce/r e , 
the true interaction can only be approximated by con- 
tact interactions up to order ("fl) 2 . If a higher accuracy 
is required, the van der Waals tail has to be taken into 
account explicitly [lj. For positive effective range, we 
pcrfomed all calculations with a cutoff well above the 
breakdown scale of the theory A 3> 1/V e - F° r negative 
effective range the calculation at next-to-leading order 
(NLO) has to be performed with cutoffs A ~ l/r e . This 
is due to a cancellation between even and uneven orders 
in the expanded two-body propagator at values of the 



loop momentum ij > 1 /r e [14j . In uneven orders of the 
expansion in r e this would lead to an incorrect renormal- 
ization if the cutoff was taken much larger than l/r e . In 
next-to-next-to- leading order (NNLO), the term propor- 
tional to (7?" e ) 2 dominates the expansion of the two-body 
propagator for large loop momenta and the problem is 
not present. 

We calculate the effective range corrections to the cor- 
relation between a and a a d two different ways: (i) in the 
subtractive renormalization scheme of Ref. [14j where the 
three-body force term Ga(A) does not explicitly appear 
and (ii) by performing a calculation at fixed cutoff A and 
varying G3. 2 Both methods agree to very high precision. 



C. Correlation between a and a a d 

The results for the correlation between a and a a d for 
positive effective range r e > are shown in Fig. [1] for 
two different values of the effective range. The solid 
line gives the leading order result, while the dashed and 
dash-dotted lines give the NLO and NNLO results, re- 
spectively. The light (dark) curves are for "fr e = 0.1 
(jr e = 0.3). For large absolute values of a ac t and 
7r e = 0.1 the shift from LO to NLO is accidentally 
very small. However, the shift from NLO to NNLO is 
of the expected size of (7r e ) 2 ~ 0.01 and is larger than 
from LO to NLO. A similar observation holds for the case 




2 Since proper renormalization has been explicitly verified, practi- 
cal calculations can be carried out at fixed cutoff A. 



FIG. 1: Correlation between a and a a d for r e > 0. The solid 
line gives the leading order result, while the dashed (dash- 
dotted) lines give NLO (NNLO). The light (dark) curves are 
for "/r e —0.1 (~/r e — 0.3). 



7r e = 0.3. We observe that the shift from LO to NLO 
varies strongly in size and is sometimes smaller than ex- 
pected by naive dimensional analysis. However, the shift 
from NLO to NNLO is generally of the expected order of 
magnitude. The smallness of the NLO corrections can be 
understood as a cancellation between two different con- 
tributions to this correction. The two contributions are 
proportional to jr e and «r e , respectively, where k is the 
typical momentum scale of the process. Furthermore, 
it is known from experience [l[ that at LO observables 
are often described better than expected from the power 
counting once the exact pole position of the two-body 
propagator is reproduced. As a consequence, the shifts 
in observables from LO to NLO can be very small and of a 
size comparable to the corresponding shifts from NLO to 
NNLO. The minimum in the correlation curve is moved 
to larger values of a a d- The size of this shift is of the 
order 7r e from LO to NLO, and of the order (7r e ) 2 from 
NLO to NNLO. 

We emphasize that an error estimate determined by 
the truncation error of the EFT expansion can be given 
from dimensional analysis within the EFT approach. 
However, this does not necessarily imply that the shift 
from one order to the next has to be exactly equal to 
this error estimate. Firstly, the error estimate is only 
accurate up to a number of 0(1) which depends on the 
details of the problem and cannot be estimated a priori. 
Secondly, if there are competing dynamical contributions 
at a given order, cancellations may lead to significantly 
smaller corrections as can be seen in our results. 

The results for the correlation between a and a a d for 
negative effective range r e < are shown in Fig. [2] for 
two different values of the effective range. Again, the 
solid line gives the leading order result, while the dashed 
and dash-dotted lines give the NLO and NNLO results, 
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FIG. 2: Correlation between a and a a d for r e < 0. The solid 
line gives the leading order result, while the dashed (dash- 
dotted) lines give NLO (NNLO). The light (dark) curves are 
for "/r e = —0.1 ("/r e — —0.3). 



respectively. The light (dark) curves are for jr e = —0.1 
(l r e — —0.3). For jr e — —0.1, we find that the cor- 
rections scale as expected. In the case jr e = —0.3, the 
behavior at NNLO resembles the case for positive effec- 
tive range. The minimum is shifted to smaller values 
of a a d- From LO to NLO the size of the corresponding 
shift is of the order jr e and from NLO to NNLO it is 
approximately of the order (7r e ) 2 . 

For negative effective range, the maximum value of 
amj 4 /% in NLO is very sensitive to the value of r e , while 
at NNLO this sensitivity disappears and the maximum 
value is close to the LO result. For r e > 0, the sensitivity 
persists at NNLO. It is also evident that the convergence 
of the EFT expansion may be slow for large absolute 
values of a a d- For smaller \a a d\, rapid convergence can 
be achieved for larger values of jr e . Indeed it will be 
shown below that for 4 He atoms the encountered values 
of "fr e are well within the range of applicability of the 
EFT expansion. The correlation curve has a pronounced 
minimum and is not completely symmetric around this 
minimum value which lies in the interval 1 < 7a™ 11 < 2. 
As a consequence, there will always be 2 solutions for a a d 
if a is used as three-body input. 



D. Application to 4 He atoms 

As an illustration, we apply our results to 4 He atoms 
for which the condition jl -C 1 is well satisfied [l], UH, [l(| . 
The interatomic potential between two 4 He atoms does 
not support any deep dimers, such that our calculations 
are applicable. The binding energies of the 4 He trimers 
have been calculated accurately for a number of dif- 
ferent model potentials for the interaction between two 
4 He atoms. (See, e.g., Refs. pjl EH and references 



therein.) The 4 He dimer and trim er g round states have 
been observed in experiment [20l . |2l| . but the excited 
state of the trimer has not yet been seen. Due to the 
limited experimental information in the three-body sec- 
tor, we apply our effective theory to the TTY [22| and 
HFD-B3-FCI1 [25[ potentials. The input parameters for 
the effective theory are the dimer binding energy and 
the effective range in the two-body sector plus one three- 
body observable. Once these parameters are specified, all 
other three-body observables can be predicted to NNLO 
in the EFT expansion. In the following, we use the ef- 
fective theory to predict a if a a d is known for a given 
potential and vice versa. 

We consider first the TTY potential [II]. The scat- 
tering length for the TTY potential is a = 100.01 A. 
This is much larger than its effective range r e = 7.329 A, 
which is comparable to the van der Waals length scale 
Kdw = 5-398 A [H The binding energy of the 4 He 
dimer for the TTY potential is B 2 = 1.30962 mK which 
is small compared to the natural low-energy scale for 4 He 
atoms: £vdw ~ 420 mK. This energy B 2 corresponds to 
a pole position 7 = 0.01040 A -1 . The atom-dimer scat- 
tering length for the TTY potential was calculated in 
Ref. [23 as a ad = (1.205 ± 0.001)7" 1 . Using this infor- 
mation as the three-body input and the effective range 
for the TTY potential, r e j — 0.076, we can obtain the 
recombination constant a which has not been calculated 
for the TTY potential. Our results are given in Table HI 



TABLE I: Recombination constant a in units of h/(m , y 4 ) for 
the TTY potential [2^]. The errors given correspond to the 
uncertainty in the input quantity a a d- (Note that % 2 /m = 
12.1194 KA 2 for 4 He atoms.) 





LO 


NLO 


NNLO 


a [ft/(m 7 4 )] 


2.791 ±0.013 


3.899 ±0.015 


3.778 ± 0.014 



The values are numerically accurate to the number of 
digits given. Even at NNLO the numerical errors in the 
input parameters are typically negligible compared to the 
truncation error of the effective theory. As an illustra- 
tion, we have propagated the numerical error in the in- 
put parameter a a d to the results in Table HI It should be 
emphasized that the physical accuracy is determined by 
the truncation error of the effective theory. The effective 
theory error estimates are 10%, 1%, and 0.1% times a 
number of 0(1) at LO, NLO and NNLO, respectively. 
The actual difference between LO and NLO and NLO 
and NNLO in our results is 40% and 4%, respectively, 
but still in agreement with the above expectation. From 
this pattern we expect our NNLO result to be accurate to 
0.4%. This slightly larger error than expected is due to 
the value of amj 4 /Ti close to the minimum in the corre- 
lation curve, where small corrections are relatively more 
important. In Ref. (l5j . a was calculated at LO using the 
trimer excited state energy instead of a a d as the three- 
body input and the value a — 2.9h/(mj' 1 ) was found. 
This result is in good agreement with our value at LO. 
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For the HFD-B3-FCI1 potential [Hj], the situation is 
the opposite and a ac i is unknown. However, the recom- 
bination constant a = 12 • 10~ 29 cm 6 /s has been cal- 
culated [26[ . The scattering length for this potential is 
a = 91.0 A, which is again much larger than the effective 
range r e = 7.291 A and the van der Waals length scale 
4dw = 5.398 A [l|,[23j]. The position of the dimer pole is 
7 = 0.01149 A^ 1 . From this information, we can extract 
r e 7 = 0.084 and amj 4 /h = 1.32. Because of the shape 
of the curve in Fig. [TJ we have now two solutions for a a d 
as shown in Table HIl 

TABLE II: The two solutions for the atom-dimer scattering 
length a ad in units of 7 _1 for the HFD-B3-FCI1 potential 
(Note that h 2 /m = 12.1194 KA 2 for 4 He atoms.) 





LO 


NLO 


NNLO 




(solution 1) 


1.341 


1.434 


1.424 


aadh' 1 } 


(solution 2) 


1.931 


2.034 


2.024 



Using the information on a ac i from the TTY poten- 
tial above, we can identify solution 1 as the appropriate 
one for 4 He atoms. The effective theory error estimates 
are 10%, 1%, and 0.1% times a number of 0(1) at LO, 
NLO, and NNLO, respectively. In this case, the differ- 
ence between LO and NLO as compared to the difference 
between NLO and NNLO is a factor of 0.7 smaller than 
the above estimate and thus in good agreement with the 
expectation. From this pattern we expect the error of our 
NNLO result to be 0.1 %. In Ref. [l[, a a d was calculated 
at LO using the trimer excited state energy instead of a 
as the three-body input and the value -fa a d = 1.4 was 
found. Again, this result is in good agreement with our 
value at LO. The deviation between results using differ- 
ent three-body inputs gives a simple estimate of higher 
order corrections. For a a d, this estimate works well. In 
the calculation of a discussed above, however, this esti- 
mate is misleading because the value of a is close to the 
minimum of the correlation curve. 



E. Conclusions 

In this work we have computed the effective range cor- 
rection to the recombination coefficient a to second order 
in the expansion in jl. Up to this order, only the two- 
body scattering length a and effective range r e plus one 
three-body parameter are required as input for the ef- 



fective theory [14J. We have calculated the correlation 
between the atom-dimer scattering length a a d and the 
recombination coefficient a at LO, NLO and NNLO, and 
studied the convergence of the effective theory expan- 
sion in 7L We found that the convergence of the EFT 
expansion may be slow for large absolute values of a a d- 
For natural values |7(w| ^ 3, rapid convergence could 
be achieved for "/r e < 0.3. In the same manner, the 
correlation between a and other low-energy three-body 
observables, such as Efimov trimer energies, could be cal- 
culated. As a test, we have applied our results to 4 Hc 
atoms. We have used the TTY H and HFD-B3-FCI1 
potentials [25[ as input and predicted, using the effective 
theory, three-body observables that have not yet been 
calculated for these potentials. The convergence of the 
expansion in 7? for He atoms was found to be in agree- 
ment with the a priori expectation. 

Further efforts should be devoted to the inclusion of 
effective range corrections for recombination processes in 
ultracold gases of alkali atoms. Such atoms have many 
deep two-body bound states which modify the recom- 
bination into the shallow dimer and provide additional 
channels for recombination into the deep states. At LO in 
the expansion in 7/ these effects have already been calcu- 
lated in this effective theory [l], 0, H, [l(| • An extension to 
NNLO and the inclusion of temperature dependence (see 
Refs. [I3,[28|,[2!|) should allow for a precise description of 
experimental loss rates such as obtained by the Innsbruck 
group for 133 Cs atoms [30| . The recombination of cold he- 
lium atoms at nonzero collision energies has already been 
considered by Suno and coworkers [31] . A further step 
to increase the precision with which these experiments 
can be described theoretically is to include the recom- 
bination into Efimov trimers and dimer-dimer scattering 
resonances [H, [H| . Recent results in the four-body sec- 
tor give hope that this goal can be achieved within the 
present effective theory framework [34], Hf| HH • 
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